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ABSTRACT 



o 

Aims. The goal of this work is to infer the star formation properties and the mass assembly process of high redshift (0.3 < z < 2.5) 
^ galaxies from their IR emission using the 24 yum band of MlPS-Spitzer 

Methods. We used an updated version of the GOODS-MUSIC catalog, which has multiwavelength coverage from 0.3 to 24 ^m and 
I ' either spectroscopic or accurate photometric redshifts. We describe how the catalog has been extended by the addition of mid-IR fluxes 

O ' derived from the MIPS 24 yum image. We compared two different estimators of the star formation rate (SFR hereafter). One is the total 

infrared emission derived from 24 /jm, estimated using both synthetic and empirical IR templates. The other one is a multiwavelength 
C/2 , fit to the full galaxy SED, which automatically accounts for dust reddening and age-star formation activity degeneracies. For both 

estimates, we computed the SFR density and the specific SFR. 

Results. We show that the two SFR indicators are roughly consistent, once the uncertainties involved are taken into account. However, 
they show a systematic trend, IR-based estimates exceeding the fit-based ones as the star formation rate increases. With this new 
catalog, we show that: a) at z > 0.3, the star formation rate is correlated well with stellar mass, and this relationship seems to steepen 
with redshift if one relies on IR-based estimates of the SFR; h) the contribution to the global SFRD by massive galaxies increases 
with redshift up to =^ 2.5, more rapidly than for galaxies of lower mass, but appears to flatten at higher z; c) despite this increase, the 
most important contributors to the SFRD at any z are galaxies of about, or immediately lower than, the characteristic stellar mass; 
d) at z - 2, massive galaxies are actively star-forming, with a median SFR ^ 300MQyr"'. During this epoch, our targeted galaxies 
assemble a substantial part of their final stellar mass; e) the specific SFR (SSFR) shows a clear bimodal distribution. 
lO ' Conclusions. The analysis of the SFR density and the SSFR seems to support the downsizing scenario, according to which high mass 

' galaxies have formed their stars earlier and more rapidly than their low mass counterparts. A comparison with renditions of theoretical 

0^ . simulations of galaxy formation and evolution indicates that these models follow the global increase in the SSFR with redshift and 

. predict the existence of quiescent galaxies even at z > 1.5. However, the average SSFR is systematically underpredicted by all models 

considered. 



Key words. Galaxies: evolution - Galaxies: high-redshift - Galaxies: fundamental parameters - Galaxies: photometry - Galaxies: 
starburst 

1 . Introduction ing to which massive galaxies evolve rapidly up to z ~ 1.5 and 

then more gradually until the present epoch, while less massive 
Answering the basic questions about the birth, formation, mass galaxies continue to evolve, implies that massive galaxi es must 
build-up, and evolution of galaxies throughout cosmic time are jjave already formed by z ~ 1 .5. Several groups (e.g., Fab eTltd] 
some of the major goals of observational extragalactic astron- ^5o||Brown et ^112007 ) smdied the evolution of massive galax- 

ies at z < 1, and their migration from the blue cloud to the red 
In the past few years, this issue has been approached sequence. We no te that optical observations (iBell et allllOOl 
with two different methods. Many previous works measured IZucca et alJl2006h suggest that the number of massive galaxies, 
a rapid evolution in the stellar mass density between z ~ as well as the stellar mass on the red sequence, has nearly dou- 
1 and z ~ 3 (jPickinson et al. 2003; Fontana et al. 2003', bled since z ~ 1, in qualitative agreement with the hierarchical 
2004; Gl azebroo k et al. 2004; Drory et al. 2004; Fontana et si. merging scenario 
2006; Rudnick et al. 2006; Papo vich et al.l2006l;iYan et al.l2006[ 

Pozzetti et al...200 7) and demonstrated that a substantial fraction A parallel line of study has analyzed the rate at which galax- 



(30-50%) of the stellar mass formed during this epoch. The dif- ies form stars during different epochs and shown that they ex- 

ferential evolution in the galaxy stellar mass function, accord- perie nce an ex tremely active phase in the same redshift range 

(e.g .. iLillvetal . 1996; Madauetal. 1996; Steiddet al. 19^ 

Send offprint requests to: R Santini, e-mail: Ho pkinsI |2004 iHopkins & Beaco m 2006; iDaddi et alj|2007t 



santiniOoa-roma. inaf . it Galaxies appear to form their stars following the so-called down- 
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sizing scenario, in which the star formation shifts from high mass 
to low mass gala xies as re dshift decreases. This picture was 
first introduced by ICowieT t al. (1996), who studied the evolu- 
tion in the Ks band luminosity function with redshift, and pro- 
poses that the most massive galaxies assemble their mass both 
earlier and more quickly than their lower mass counterparts, 
which, in contrast, continue to form stars until recent epochs. 



Later on, many other groups (e 
Fontana et al. 2003; Feulner et al 



Brinchmann & Ellis 20QQ: 



20051; Perez-Gonzalez et al] 



2005; iPapovich et al.,,2006: .Damen et al.ll2009l) derived confir 



mations of a downsizing behaviour from the study of the specific 
star formation rate, defined to be the star formation rate per unit 
mass, at diff erent redshifts. Ho wever, we note that deep radio 
observations dPunne et al.ll2009l) appear to conflict with this sce- 
nario. The downsizing picture also appears to contradict the hi- 
erarchical growth scenario in which the most massive structures 
that we see today are produced by merging processes between 
smaller structures inside large-scale overdensities and collapsed 
when the Universe was far younger than today. 

To re produce this early formation of massive galaxies (see 
iThomas et al. 2005) that are akeady "red and dead" at high 
z, theoretical models had to introduce very efficient processes 
of star formation and its suppression by means of active 
galactic nuclei and supernovae quenching of c ooling flows 
(Menci et al.' 2006; Kitzbichler & White 2007; Bower et al.' 
2006; Croton et al. 2006; Nagamine etal. 2006; Monaco et aL 
20071). gravita t ional heating dKhochfar & Ostriker 20081 
Johansson et al.l l2009l) . or shock heating dPekel & Birnboir 



20061) . These models differ slightly in their predictions mainly 



because they adopt different processes to shut down the star 
formation. 

Both stellar masses and star formation rate estimates are af- 
fected by a number of uncertainties. The measure of the star 
formation rate (SFR) is especially difficult to handle. The high 
amount of energy produced by newly born stars is emitted 
throughout the galaxy spectral energy distribution (SED), from 
X-rays to radio frequencies. Since we are of course unable to 
directly measure the total light emitted by young and massive 
stars, calibration factors and corrections a re apphed to estimate 
its va l ue for any of th ese frequency ranges (lKennicuttll998l:lBeil 
120031; ICalzettilf2008l) . One of the most commonly used estima- 
tors is the UV rest-frame band, where young and massive stars 
emit most of their light. However, dust absorbs, reprocesses, and 
re-radiates UV photons at near-to-far IR wavelengths. Hence, 
the reliability of UV luminosity as a SFR tracer depends on 
large and uncertain corrections rel ying upon the dust properties, 
which are n ot yet clearly known dCalzetti et al.l 119941; ICalzettil 
119971 12 001|). Moreover, the UV-upturn at A shortward of 2500 
A (e.g.. iHarTet al . 2007). especially in elliptical galaxies, can 
potentially bias the SFR estimate at very low redshift. Since 
the most intense star formation episodes are expected to occur 
in dusty regions, most of the power originating in star-forming 
(SF) galaxies is emitted in this wavelength range, and the dust 
emission peak is the dominant component of SF galaxies SEDs 
dAdelberger & Steidelll2000HCalzetti et al.ll2000h . Thus, a popu- 
lar approach consists of adopting a conversion between the to- 
tal emitted IR luminosity (Lir hereafter) and a star f ormation 
rate e stimation that is unaffected by dust obscuration dKennicutti 
II998I) . 

The total infrared luminosity is generally estimated by com- 
paring observed SEDs and synthetic templates, although empiri- 
cal conversions have s ometimes been used (Takeuchi et al. 2005; 
iBavouzet et al.l l2008l). A variety of different hbraries are used 
for fliis purpose (e.g.. lChary & Elbazll2001l;lDale & Heloul2002l; 



iLagache et al.ll2003t ISiebenmorgen & Krugelll2007l and so on). 
A notable problem for the reliability of IR-based SFR tracers 
concerns obscured AGNs. In these objects the IR emission is 
generated by matter accretion onto a central black hole rather 
than dust heating by young stars. 

In this paper, we use the GOODS-MUSIC catalog to inves- 
tigate properties of star-forming galaxies up to redshift 2.5 and 
infer the mass assembly process from thek mid-IR emission. The 
paper is organized as follows. In Sect. 2, we recall the basic fea- 
tures of the GOODS-MUSIC dataset and explain the innovations 
concerning its latest version, and we explain how it has been 
updated with the addition of the 24 jjm photometric band. We 
derive and compare star formation rates from IR- and fit-based 
estimators in Sect. 3. In Sect. 4, we present a study of the mass 
assembly process in the high redshift Universe and a compar- 
ison with theoretical model predictions. We finally summarize 
our work and our conclusions in Sect. 5. In Appendix A, we de- 
scribe in more detail how we convert mid-IR fluxes into total 
infrared luminosities and compare the different templates used, 
and in Appendix B we present the error analysis performed on 
the fit-based SFR estimates. 

Throughout t his work, unless stated otherwise, we assume 
a ISalpeterl dl955h initial mass function (IMF) and adopt the A- 
CDM concordance cosmological model (Hq - 70 km/s/Mpc, Q.m 
= 0.3 and Qa = 0.7). 



2. The data sample 

2.1. The new GOODS-MUSIC sample 

We present and use an updated version of the multicolour 
GOODS-MUSIC sample (GOODS Multicolour Southern 
Infrared Catalog; [Crazian et al.ll2006l) . extracted from the pub - 
Uc data of the GOODS-South survey dGiavahsco et al.ll2004 . 
In the following, we shall refer to this version of the catalog as 
GOODS-MUSIC v2, to differentiate it from the former public 
version, which is named vl hereafter. The new version is also 
made publicly available Q. 

The 15-bands multiwavelength coverage ranges from 0.35 to 
24 fim, as a result of the combination of images from different 
instruments (2.2ESO, VLT-VIMOS, ACS-HST, VLT-ISAAC, 
Spitzer-IRAC, Spitzer-MIPS). The catalog covers an area of ~ 
143.2 arcmin^ located in the Chandra Deep Field South and con- 
sists of 15 208 sources. After culling Galactic stars, it contains 
14999 objects selected in either the z band or the band or at 
4.5 fim. 

The whole catalog has been cross-correlated with spectro- 
scopic catalogs available to date, and a spectroscopic redshift 
has been assigned to ~12 % of all sources. For all other objects, 
we have computed well-calibrated photometric redshifts using a 
standard minimization technique for a large set of synthetic 
spectral templates. 

The previous version of the catalog and procedures adopted 
to determine the photometric redshifts an d physical properties 
of each object were described at length in iGrazian et al.l (I2006h 
and Fontana et al. (2006). With respect to the previous catalog, 
we have performed a set of improvements to the optical-near-IR 
data, the major ones being: 



' The catalog is available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ It is also possible to 
download the catalog at the WEB site http://lbc . mporzio . astro .it/goods , 
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- In addition to objects selected in the ACS z and in the ISAAC 
Ks bands, we have also included objects selected from the 
IRAC 4.5 jum image, hence including sources detected at 4.5 
//m but very faint or undetected even in Ks band. A full de- 
scription of these objects is beyond the scope of the present 
paper and will be presented elsewhere. 

- We have revised the photometry in the four IRAC bands us- 
ing an updated version of the PSF-matching kernels, as re- 
leased by the Spitzer Science Center. As a consequence, we 
also adopted a larger kernel, to fully account for the large 
tails of the IRAC PSFs. 

- We have adopted a revised procedure for estimating the 
background in the IRAC images. Using the objects' posi- 
tions and IRAC fluxes from the GOODS-MUSIC vl catalog, 
we have created realistic simulated images in the four IRAC 
bands by smoothing sources to the nominal IRAC PSFs. An 
accurate background estimation has been performed by sub- 
tracting these simulated images from the original ones and 
by linearly interpolating the residual emission. Since the av- 
erage value of the apparent IRAC background is negative, 
this has led to an increase in the adopted background, with 
respect to the GOODS-MUSIC vl version. 

- These two changes have modified the IRAC photometry. 
Because of the new kernels, the brightest objects have a 
higher flux and a typical offset in magnitude of 0.23, 0.14, 
0.22, 0.35, respectively in the 3.6, 4.5, 5.8, and 8 yum im- 
ages, with respect to the G OODS-MUSIC v l catalog, which 
agrees with the analysis of lWuvts et alj (1200 8). This effect is 
largely mitigated for fainter galaxies, since the higher back- 
ground now adopted has led to an underestimate of their 
fluxes. 

- Overall, the revised IRAC photometry has a modest impact 
on the estimate of photometric redshifts, since {Zphm v2 - 
Zphot yi) ~ 0.01 + 0.16 for the whole sample and ~ 0.01 + 
0.03 when a 3cr-clipping analysis is performed (see also 
IWuvts etal] ( l2008) ). 

- A more informative test of the photometric redshift accu- 
racy comes from the enlargement of the sample of galax- 
ies with spectroscopic redshifts, which w e obtained by 
adding new spectra from public surveys dVanzella et alj 
l2008tlPopesso et al.1 12009) . In addition, we have also had 
access to the spectra of the GMASS survey dCimatti et al.l 
|2008|) . prior to their publication. The final sample now in- 
cludes 1888 galaxies, three ti mes larger than the spectro- 
scopic sample in iGrazian et al.i (l2006l) . The additional spec- 
tra are mostly relative to galaxies that are both fainter and 
at higher redshift than in the original sample. Without sig- 
nificant refinements in the adopted templates, we then find 
that the absolute scatter |Az| - \zspe - Zphotl/i^ + Zspe) has a 
slightly larger average value. Quantitatively, the average ab- 
solute scatter is now (|Az|) = 0.06, instead of 0.045 obtained 
for the GOODS-MUSIC vl catalog. However, when only the 
brightest galaxies are considered, we find comparable values 
with respect to Grazian et al. (2006) «|Az|) = 0.043). We 
have verified that this is due to an increased number of out- 
liers, as shown by a 3cr-clipping analysis, which provides 
(|Az|) - 0.027 and 0.032 for the complete datasets of vl and 
v2 catalogs, respectively. 

- We have removed Galactic stars and performed a more care- 
ful selection of the galaxy sample to identify AGN sources. 
For the latter, we have first removed all objects whose spec- 
tra show AGN features. Then, w e have cross-correl ated our 
catalog with the X-ray catalog of Brusa et al] d2009'). and re- 
moved all X-ray detected sources whose flux is dominated 



by an unresolved central source. These sources typically 
have spectra classified as narrow-line AGNs. The optical 
morphologies of all remaining X-ray sources do not show 
a dominating central point like source, and - where avail- 
able - have typically spectra classified as emission-line star- 
forming galaxies. These objects have been retained in our 
galaxy sample. 

The major new ingredient of this new version of the 
GOODS-MUSIC catalog, however, is the inclusion of the 
24 fim photometry for all galaxies in the sample, which is the 
main focus of the present paper We describe the adopted proce- 
dures and results in the following. 

2.2. MIPS 24 fim catalog 

We have extended the GOODS-MUSIC catalog with the addi- 
tion of the mid-IR fluxes derived from the public 24 //m im- 
age of the Multiband Imager Photometer for Spitzer (MIPS, 
Rie ke et ani20 04) onboard the Spitzer Space Telescope. As for 
the IRAC images, given the very large PSF of this image (~ 5.2 
arcsec), to properly detect and de-blend objects we had to em- 
ploy a PSF-mat ching techniq ue, which is performed by the soft- 
ware ConvPhot (iDe Santis et al. 2007). This algorithm measures 
colours from two images of different qualities by exploiting the 
spatial and morphological information contained in the higher 
resolution image. 

When applying ConvPhot to our case, each object was ex- 
tracted from the high resolution z band ACS-HST (PSF ~ 0.12") 
image, which was used as a prior to extract the objects' posi- 
tions, filtered with a convolution kernel, and finally scaled by a 

minimization over all the image pixels to match the intensity 
in the MIPS image. To fully use the positional information of the 
ACS images, and maintain consistency with the z-selected sam- 
ple, we used the z band ACS-HST data as a prior, augmented 
by artificial objects placed where galaxies detected only in Kg 
or 4.5 jum were located. In our case, the MlPS-Spitzer and ACS- 
HST images have pixel scales of 1 .2 and 0.03 arcsec/pixel, which 
means that it is impracticable to use ConvPhot even with fast 
workstations. To make the computation feasible, we rebinned 
the ACS detection image by a factor 8x8 (0.24 arcsec/pixel). In 
regions where the crowding of the z-detected sources was signifi- 
cant, the fit may become unconstrained or degenerate because of 
the large size of the MIPS PSF. To prevent this, we placed an ad- 
ditional constraint on the fitted fluxes that must be non-negative. 
For the objects whose flux is forced to be zero, we provide an up- 
per limit derived from the analysis of the mean rms in the object 
area. 

The behaviour of ConvPhot in these extreme applications 
was tested with sever al simulations, which were described by 
iDe Santis et al. I d2007h . These tests indicated that the estimated 
magnitudes are biased by neither the different qualities of the 
two images nor the undersampling of the high resolution image 
after rebinning. 

Finally, from visual inspection of sources with unusual 
colours, we removed ~ 30 objects from the catalog whose flux 
was incorrectly assigned. By examination of the residuals, we 
also verified that there is no significantly bright source in the 
MIPS image apart from the ones considered. We detected 3313 
(~ 22% of total) detected objects and 1 1 841 (~ 78%) Icr upper 
limits. 

Despite the validation tests providing satisfying results, we 
emphasize that the intrinsic limitations due to the poor resolution 
of the MIPS image cannot be completely overcome, in particu- 
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lar for sources blended even in the detection image, i.e., sources 
whose profiles overlap in the ACS z image. In this case, the sep- 
aration between the two - or more - objects can become even 
smaller than the positional accuracy of the MIPS images, and the 
association between the z-detected and the MIPS sources relies 
on the accuracy of the astrometric solution. To check for possi- 
ble misidentifications, we associated each source with a flag to 
quantify the number of possible contaminants, at different dis- 
tances. This flag was attached to the public catalog and we cau- 
tion the user to check for possible systematic errors. The results 
that we shall present in the following are insensitive to the inclu- 
sion of the most blended sources, which have therefore not been 
removed. 

It may be interesting to observe the effect of this proce- 
dure on the 24 yum number counts, shown in Fig. [T] We present 
both the counts derived by ConvPhot as well as those obtained 
by a SExtractor catalog, which was produced by adopting cor- 
rected aperture magnitudes at 6 arcsec. The blue shaded re- 
gion shows the 24 yum counts from lPapovich et a l. (2004). They 
counted sources in five different fields (CDFS among them) us- 
ing DAOPHOT software. The green solid line (err or bars are not 
large enough to be seen) represents IShupe et al.l (|2b08) counts 
in the SWIRE field performed by SExtractor. ConvPhot (black 
dots) provides consistent results with both the previous work of 
other authors to the lowest fluxes and our SExtractor catalog (red 
circles) at F > 100 fiJy. 

The agreement at bright fluxes confirms th e results of the 
simulations and validation tests presented in iDe Santis et al.l 
(1200 7) and shows that the fluxes estimated by ConvPhot 
agree with those estimated by an independent detection with 
SExtractor for sources that are not severely blended. For the 
blended fraction of the objects, however, the fluxes estimated by 
ConvPhot can be slightly lower, since part of the flux is identi- 
fied with the fainter contaminants, which are not detected in the 
24 jum image alone. 

It is more interesting to consider the behaviour at faint fluxes, 
where the a priori knowledge of the object position, because of 
the use of the z, and 4.5 ;um images for the detection, al- 
lows us to obtain flux estimate at much fainter limits, reduc- 
ing the effects of blending and confusion. As expected, indeed, 
SExtractor counts decline at ~ 100 fxiy, where the confusion limit 
prevents the detection of fainter sources, while ConvPhot allows 
us to complete our detection to even deeper limits. ConvPhot 
number counts present a double slope, with a break point lo- 
cated at ~100 yuJy, which we consider to be an intrinsic prop- 
erty of the sample. The slope and the normalization at the faint 
end agree with the estimates of lPapovich et al.l (|2004|) . who care- 
fully computed a correction for the incompleteness due to poor 
resolution at faint limits. Papovich et al. (2004) and subsequent 
papers (e.g., [L e Floc'h et al. 2005; Perez-Gonzalez et al. 2005t 
iPapovich et"al..2006;,Marcillac et al..2006;,Bell et al.,200 7) esti- 
mated that the source detection in MIPS-CDFS is 80% complete 
at 83 yuJy or so, which typically corresponds to S IN ~ 2% - 30 
in our fitting procedure. For consistency with these works, we 
therefore distinguish between objects with fluxes above this limit 
and those detected at lower S IN, to fluxes as low as ~20 Jy, 
which corresponds to our flux count limit. The median S IN at 
this flux limit is ~ 6-7, although a tail is present at lower values 
of S IN that is caused by source blending. Nevertheless, this tail 
includes only a small number of objects. 



T3 
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Fig. 1. 24 yum flux number counts measured by ConvPhot (black 
filled circles) and SExtractor (red empty circles). Error bars have 
been computed as the square root of the number of objec ts in 
each bin. We compared our results tO iPapovich et al.l ( |2004|) and 
IShupe et al] ([2008) 24 fim counts (blue shaded region and green 
line, respectively). 



one released by the GOODS Teanfl dCharvl |2007|) . who 
adopted a similar source extraction technique. The over- 
all agreement was good with a small offset in both 

cases ((-F24//m goods -music IP lAiim FIREWORKS 

) ~ 1.2 and 

{FlAiim GOODS-MUSIC I FlAiim DRi) ~ 0.91). 

2.3. The data selection 

As pointed out above, we removed from our catalog Galactic 
stars and both spectroscopic or X-ray detected AGN sources. 
Moreover, we only consider the redshift range 0.3 - 2.5. From 
this sample, we consider the following two subsamples: 

- the purely /fj-selected sample, subsample A, consists of 
2602 galaxies. Of these, 981 (376) are 24 jum detections with 
F24/jm > 20 (83) fiiy and 1621 have been assigned an upper 
limit; 

- subsample B was created by performing the following cuts: 
z <26otKs < 23.5 or m4 5 < 23.2. It includes 7909 galaxies, 
of which 1165 (413) are 24 fim detections with F24^,„ > 20 
(83) //Jy and 6744 are upper limits. 



3. Comparison between SFR indicators 

In this section, we present different methods to estimate the star 
formation rate and compare the results obtained from the differ- 
ent indicators. 



Finally, we cross-correlated our 24 ^m cata log with 2 GOODS-South MIPS 24 micron source list v0.91 from GOODS 
both the FIREWORKS catalog (IWuvts et al.1 12008|) and the data release (DR3). 



p. Sandni et al.: Star formation and mass assembly in high redshift galaxies 



5 



3.1. SFR estimators: IR-, fit- and UV-based SFR 

Based on the assumption that most of the photons originating 
in newly formed stars are absorbed and re-emitted by dust, the 
mid-IR emission is in principle the most sensitive tracer of the 
star formation rate. In addition, a small fraction of unabsorbed 
photons will be detected at UV wavelengths. A widely used SFR 
indicator is therefore based on a combination of IR and UV lu- 
minosity, which supplies complementary knowledge of the star 
formation process (Iglesias-Paramo et al. 2006; Calz etti et al.l 
120071) . For 24 yum detected sources, we estimated the instanta- 
neo us SFR us i ng the same calibration as iPapovich et al.l (l2007l) 
and lBelletall(l2005h : 



SFR 



IR+UV 



/Mayr 



1.8x 10""'xLto,/Lc 



(2.2 X Luv + Lir) 



(1) 



(2) 



We computed Lir by fitting 24 /im emission to lDale & Heloul 
( I2OO2I) (DH hereafter) synthetic templates, which are widely 
adopted in the literature. In Appendix A, we compare these val- 
ues wit h the corresponding L i r predicted by differ ent model li- 
braries (ICharv&Elbazll2001l; iPoUetta et aLll2007h . We include 
the rest-frame UV luminosity, uncorrected for extinction, de- 
rived from the SED fitting technique, Luv = 1.5 x L2700A; al- 
though often negligible, this can account for the contribution 

from young u nobscured stars. 

Following lPapovich et al.l(l2007l) . we then applied a lowering 
correction to the estimate obtained from Eq. [1] They found that 
the 24 yum flux, fitted with the same DH library, overestimates 
the SFR with respect to the case where longer wavelengths (70 
and 160 jum MIPS bands) are considered as well, and they cor- 
rected the trend using an empirical second-order polynomial. In 
Appendix A we show further confirmation of the need to ap- 
ply this correction: for bright sources, Lir estimated by syn- 
thetic models has values of up to a fact or of 10 higher tha n 



Lir predicted by the empirical libra ry of PoUetta et al. J 20071) 



■y 

Similar r esults were also published bv lBavouzet et al] (1200 8') and 
iRieke et al. (2009). In the following, we refer to the estimate in 
Eq.[T]as SFRiR+uv- 

A complementary approach to estimating the star formation 
rate, as well as other galaxy physical properties (e.g., mass, age, 
dust extinction), is the SED fitting. A grid of spectral templates 
is computed from standard spectral synthesis models, and the 
expected magnitudes in our filter set are calculated. The derived 
template library is compared with the available photometry and 
the best-fit model template was adopted according to a min- 
imization. During the fitting process, the redshift is fixed to its 
spectroscopic or photometric value. The physical parameters as- 
sociated with each galaxy are obtained from the best-fit template 
up to 5.5 fim rest-frame. This analysis assumes that the overall 
galaxy SED can be represented as a purely stellar SED, extincted 
by a single attenuation law, and that the relevant E{B-V) and ba- 
sic stellar parameters (mostly age and star formation history, but 
also metallicity) can be simultaneously recovered with a mul- 
tiwavelength fit. We note that parameter degeneracies cannot be 
completely removed, especially at high reds hift. P revious studies 
(iPapovich etal.1120011; iShaplev et al.i 120011 l2005b demonsti-ated 
that, while stellar masses are well determined, the SED fitting 
procedure does not strongly constrain star formation histories at 
high redshifts, where the uncertainties become larger due to the 
SFR-age-metallicity degeneracies. For this reason, the uncer- 
tainties associated with the SFR values estimated from the SED 
fitting are larger than those associated with the IR tracer 



In our analysis, w e estimated star formation r ates (along with 
stellar masses) using Bruzu al & Charloi ( l2003h synthetic mod- 
els, fitting the whole 14 bands of photometry (from the U band 
to 8 fim). We parameterize the star formation histories with a 
variety of exponentially declining laws (of timescales t rang- 
ing from 0.1 to 15 Gyr), metallicities (from Z - 0.02 to 
Z - 2.5 Zq) and dust extinctions (0 < E(B-V) < 1.1, with a 
Calzetti or Small Magellanic extinctio n curve). Detai l s are g iven 
in Table 1 ofLF ontana et al. (2004J), in iFontana eTal] (l2006h and 
in lGrazian et al.. (i2006t,2007.) . With respect to these recipes, the 
only difference in our method is the adoption of a minimum age 
of 0. 1 Gyr Below this value, the relation between UV luminosity 
and star formation rate changes rapidly with the age of the stellar 
population, leading to very high values of inferred SFRs. We are 
aware that exponential star formation histories may not be the 
correct choice in some cases. However, modeling in detail the 
star formation history of our galaxies is beyond the scope of the 
present paper, whose aim is to compare the star formation rates 
derived from the IR emission with those based on the widely 
used SED fitting procedures. In the following, we refer to this 
SFR estimation as SFRgt. 

We have also fitted our data using the Maraston (20051) 
(MR05) and Chariot & Bruzual (in prep., see Bruzual 2007 aibb 
(CB07), including an improved TP-AGB stars treatment. The 
stellar mass estimates inferred using these new models are pre- 
sented in Salimbeni et al. (2009) a nd are approxirn ately 0.2 dex 
lower than those computed with iBruzual & CharlotI ( 12003!) mod- 
els. However, the extrapolated SFR v alues lead to a si gnificant 
overestimate (especially in the case of lMarastonI (l2005h models) 
of the SFR density that we present in Sect. 14.21 We suspect that 
part of this discrepancy is caused by the peculiar shapes in the 
near-IR side of the spectrum, which probably caus e the worse 

tha n measured from the comparison with Bruzua ffe Charloj 
( 120031) models. Since a robust comparison among different stel- 
lar population models is beyond the scope of this paper, and 
pending further tests on the new M R05 and CB07 models, we 
continue to adopt the widely used iBruzual & Charloi (l2003h 
template library and refer to possible future work for more de- 
tails on this point. 

At z > 1.5, it is possible to obtain an independent estimate 
of the SFR using the observed Lism and the observed slope of 
the UV continuum to estimate the E(B - V), rather than the 
multiwavele ngth fit. In t he fol lowing, we use the conversions 
adopted by iDaddi et all ( |2004 . culling a sample of BzK-SF 
galaxie s from our catalog and deriving the relevant SFR us- 
ing the lDaddi et all (|2004) scaHng relations SFRuvlM^yr^ - 
Li50o/Lo, wifli Lo = 8.85 x IQ^'^erg s~^Hz^\ and A1500 = 
10 X E(B - V), where E(B - V) = 0.25(B - z + 0.1)ab- We 
refer to this SFR estimation as SFR1500. 

For clarity, SFRs lower than 0.01 M^yr^^ are assigned a 
value 0.01 Mojr ' in the following analysis. 



3.2. Comparison between IR- and fit-based tracers 

We now discuss the consistency between the different SFR esti- 
mates. In this section, we use the /Tj-selected sample (subsample 
A) to ensure a proper sampling of the full SED, e.g., to ensure 
that all the bands, or most of them, are available for the fitting 
procedure. 

We start from the redshift range 0.3 - 1.5, plotting in Fig. |2] 
the comparison between SFRir+uv and SFRfit. We include both 
the 24 yum detected and undetected galaxies, for which only up- 
per limits to SFRiR+uv can be obtained. These upper limits are 
misleading for the interpretation of the figure: since the SED 
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Fig. 2. Relation between SFRir+uv and SFRfit in three redshift 
bins for the whole /Ts-selected sample. Filled black dots are for 
the F24^m > 20 fiJy subsample; small gray dots refer to galaxies 
undetected at 24 yum, and must be regarded as upper limits. The 
blue circle identifies an obscured AGN candidate, sele cted ac- 
cordin g to a similar technique as the one presented in Fio re et aTl 
( |2008|) (see text). The red line defines the locus SFRir+uv — 
SFRfit. The typical uncertainties associated with both estimates 
are shown in the lower panel. The intrinsic parameter degen- 
eracy involved with the SED fitting procedure is responsible for 
the larger uncertainties associated with SFRgt. Units are Mgyr ' . 



fitting can reach lower nominal values for the given SFR, the 
scatter in SFRfit seems to be larger than that in SFRir+uv- 

Given the many uncertainties involved in both estimators, 
the overall consistency between them appears reassuring. Apart 
from offsets and other systematics, which we discuss below, the 
majority of galaxies are assigned a consistent SFR, and the num- 
ber of severe inconsistencies is small. These inconsistencies may 
have two different origins. 

On the one hand, galaxies with SFRfit much higher than 
SFRiR+uv can in principle be infeiTed from incoiTect fitting of 
red galaxies. For these objects, the SED fitting could erroneously 
assign a large amount of dust to an otherwise dust-free, passively 
evolving population. Despite the relatively large number of pas- 
sively evolving galaxies at z < 1 .5, the number of these misiden- 
tifications is very small at z < 1, and low even at z = 1 - 1.5. 

On the other hand, galaxies with SFRir+uv >> SFRfit can be 
obtained either when the opposite misidentification occurs (e.g., 
for dusty star-forming galaxies fitted with a passively evolv- 
ing SED) or, more interestingly, when the mid-IR emission is 



due to additional processes, not observable in the UV/optical 
regime. Typical cases are AGN emission or additional star for- 
mation activity that is completely dust enshrouded. At z < 1.5, 
these objects are again very rare in our sample. In particular, the 
SFR of the few galaxies with SFRir+uv - 2 - 6 M^yr^^ and 
SFRfit < IMojr"' at z < 1 is probably overestimated, because 
of the incorrect application of a star-forming template to a more 
quiescent galaxy (see Fig. lA.3l l. 

For a more detailed discussion of the robustness of the fit- 
based SFR estimates, the related eiTor analysis and contours, 
we refer the reader to Appendix B. 

The existing systematic trends can be appreciated by exam- 
ining Fig. [3] (a), where we plot the SFRiR+uv/SFRfit ratio for 
the 24 jum detected sample. First, we note that the scatter in the 
SFRiR+uv/SFRfit distribution widens with redshift. This is proba- 
bly due to a combination of effects: as we move to high redshifts, 
galaxies become intrinsically fainter and the rest-frame spectral 
coverage becomes narrower, making the SED fitting more un- 
certain; likewise, 24 fim based indicators could become more 
uncertain as the filter moves away from the rest-frame mid-IR 
and approaches the PAH region. 

More interestingly, a coiTelation between the 
SFRiR+uv/SFRfit ratio and SFRir+uv can be observed. If 
we focus our attention on galaxies with star formation rates of 
between 10 and 100 Mojr"', the SFRiR+uv/SFRfit distribution 
(black histograms) is centred on unity and fairly symmetric. In 
contrast, galaxies with milder activity (gray shaded histograms) 
have lower SFRiR+uv/SFRfit ratios, and we find a first hint 
of highly star-forming galaxies (red horizontally shaded his- 
togram), i.e., objects with star formation rates higher than 100 
Mojr"', that have higher SFRiR+uv/SFRfit ratios. 

The same systematic effects are seen at higher redshift. The 
upper panel of Fig. [3] (b) shows the ratio SFRiR+uv/SFRfit in 
the 1.5 - 2.5 range. It is immediately clear that the spread is 
larger than at lower redshifts, and that a large number of objects 
with SFRiR+uv/SFRfit >> 1 is observed. The observed spread 
is not surprising, since a similar disagreement, of up to two or- 
ders of magnitude in single galaxies at comparable redshi fts, has 
already been noticed by the similar analysis of Papovic ti et al.l 
(l2006l) . At these redshifts, the faintness of the galaxies and the 
large A;-corrections in the mid-IR are clearly even more effec- 
tive in increasing the noise in both estimates. The factor of ~2 
shift in the distribution of galaxies with SFRs of between 10 and 
100 Mo^r"' is consistent with the uncertainty associated with 
the spectral library used to compute Lir (see Appendix A). 

It is therefore interesting to compare these estimates with the 
pure UV-based SFR1500, described above. The relation between 
SFRiR+uv and SFRir+uv/SFRisoo rat io is shown i n the middle 
panel of Fig. [3] (b). Following D addi et al.l (|2004|) . we applied 
this recipe only to the sample of BzK-SF in the redshift range of 
interest. The scatter is significantly lower than in the upper panel, 
but the same trend for exhibiting an IR excess still appears for 
highly star-forming objects. 

To some extent, the excess in the mid-IR derived star for- 
mation rate is m ost probably caused by the presence of highly 
obscured AGNs. 'Daddi et al I ( l2007ah assumed that all the mid- 
IR excess objects drawn from the BzK-SF sample (i.e., all ob- 
jects in the lower panel of Fig.[3](b) with SFRir+uv/SFRi5oo > 3) 
are powere d by ob scured AGNs. Following a different approach, 
Fior e et al.l (l2008l) identified a population of highly obscured 
AGNs candidates in the entire Ks-selected sample, selected on 
the basis of their ver y red spect rum {FzA^iml F{R) > 1000 and 
R-K, > 4.5) (see also lDev etaiT20 08). To use the fluxes directly 
measured in our catalog, we selected the very same objects fol- 
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Fig. 3. Left panels in (a) and left top panel in (b): relation between SFRiR+uv/SFRgt ratio and SFRjr+uv in the different redshift bins 
for the F24fim > 83 //Jy (black dots) and 20 < f 24^^,,, [ju Jy] < 83 detected galaxies (gray dots); blue circles identify obscured AGN 
candidates, selected according to a similar technique as the one presented in Fiore et al. (2008) (see text). Left middle panel in (b): 
relation between SFRir+uv/SFRisoo ratio and SFRir+uv in the redshift bin 1 .5-2.5 for the subsample of BzK-SF galaxies; symbols 
are as in the other panels. Right panels in (a) and (b): black plain, gray (diagonally) shaded and red (horizontally) shaded histograms 
show SFRiR+uv/SFRfit (SFRjr+uv/ SFR1500 in the middle panel in (b)) values respectively for lO<SFRiR+uv[A/0y?""'] < 100, 
SFRiR+uv[A^03"""'] < 10 and SFRiR+uv[Moyr"^] > 100 samples. Units are M^yr'^. 



lowing the criterium F24fim/F(r) > 1000 and I-K^ > 4.5 , whic h 
we checked to be consistent with that used in lFioreetalJ(l2008l) . 
These objects are shown as blue open circles in Figs. |2] and [3] 

The presence of highly obscured AGNs could also produce 
the observed trend in the SFRiR+uv/SFRfit ratio, since they are 
expected to exist inside high mass galaxies, which are on av- 
erage highly star-forming. A correlation between the fraction 
of obscured AGNs and the stellar mass was indeed shown in 
iDaddi et all ( l2007ah . 

However, we find that the trend in the SFRiR+uv/SFRgt ra- 
tio also extends to lower star formation rates (of average values 
lower than unity) and at lower redshifts. It is therefore possible 
that it reflects a change in the intrinsic physical properties of star- 
forming galaxies. A possible explanation is related to metallicity 
effects. Galaxies with sub-sol ar metallicities have lower mid- 
IR emission (at least at 8//m, ICalzetti et al.l l2007h and higher 
UV luminosity than solar ones, for a given level of SFR. The 
observed trend can therefore be related to a metallicity trend, 
which is natural to expect given t he observed mass-me tallicity 
relation from low to high redshifts dMaiolino et alj2008l and ref- 
erences therein). Unfortunately, a direct check of the statement 
above is not feasible. Reliable metallicities cannot be inferred 
from broad-band SED fitting, and high resolution spectroscopy 
is necessary to properly distinguish between SEDs characterized 
by different lines and hence metallicities. 



Alternatively, the observed trend can be taken as evidence 
of a failure in the assumption that a single attenuation law can 
adequately model the output from a star-forming galaxy. 

4. Mass assembly and downsizing 

We now discuss the star formation properties of our galaxy sam- 
ple as a function of redshift and stellar mass. With respect to 
other surveys, our sample has the distinctive advantage of being 
selected by a multiwavelength approach. In this section, we use 
the subsample B. The Ks and m4 5 cuts ensure a proper sampling 
of highly absorbed star-forming galaxies, and hence probably 
a complete census of all galaxies with high SFR. On the other 
hand, the deep z-selected sample contains the fainter and bluer 
galaxies of both low levels of dust extinction and low star for- 
mation rate. 

In this catalog, we derive the SFR, which we refer to as IR- 
based, using the SFRir+uv estimates to all objects with F24/jm ^ 
20 fiiy, and the SFRfit to all fainter objects. We recall that the 
SFRiR+uv star formation rates are derived from the mid-IR emis- 
sion, computed with the DH s ynthetic models , and by adopt- 
ing the lowering correction of Papovich et alj ( 12007 ) for high 
SFRs. This technique was widely adopted in the Uterature, and 
we adopt it as a baseline. At the same time, we mention how 
the results would be changed by using SFRgt for all the objects. 
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Fig. 4. Left panels: relation between the star formation rate and 
the stellar mass of the GOODS-MUSIC galaxies at different red- 
shifts; large dots represent the galaxies with SFR derived from 
the 24 fim emission, while small ones are 24 yum undetected 
galaxies, with SFR derived from the SED fitting analysis; red 
solid lines show the 2cr-clipped least square fit described in the 
text; green dashed lin e in the high est redshift bin represents the 
correlation found by iDaddi et all (|2007b). Right panels: SFR 
distribution in two mass bins, 10'° - lO^M© (plain black his- 
tograms) and 10^' - IO'^Mq (shaded red histograms). 



We use the stellar masses estimates derived from the SED fitting 
analysis described above. 

Finally, to avoid bias in the IR-based SFR es timates, we 
remov e the obscured AGN candidates identified by iFiore et al] 
(I2008h from our analysis. We are aware that the removal of 
this population will cause a reduction in the star formation rate 
measured on average at high masses. AGNs a re indeed known 
to reside preferent i ally in high mass sources jBest et al.l 120051; 
iDaddi et al.|[2007al; iBrusa et alj|2009l) ; although a high fraction 
of IR emission is understood to originate from accretion pro- 
cesses, some fraction of it probably originates from star forma- 
tion. 



4. 1 . The stellar mass-SFR relation 

A direct relation between the stellar mass and the star forma- 
tion rate in high redshift galaxies was identified both at z ^ 1 
(Elbaz et al. 2007) and in a subsample of star-forming BzK-SF 
at z ^ 2 dDaddi et al.l2007bl) . To some extent, this disagrees with 
analogous trends for galaxy properties in the local Universe, 
where the m ost massive galaxies h ave very low levels of star for- 
mation rate (Heav ens et al.l l2004). We illustrate the global evo- 
lution in this relation by plotting in Fig.|4]the star formation rate 
of each galaxy as a function of its coiTesponding stellar mass. 

It is immediately clear that a trend between star formation 
rate and stellar mass is present at all redshifts. This relation is 
of course neither tight nor unique: at all redshifts there exists a 
fraction of galaxies with low star formation rates, as also shown 
by the inset histograms. The SFR distribution becomes more and 
more bimodal with decreasing redshift. Figure|4]also shows that 
at the highest redshifts only the most massive galaxies are al- 
ready experiencing a quiescent phase. We note that we are in- 
complete at low masses and SFRs because of the magnitude 
limit of the sample. Nevertheless, there is a clear trend between 
SFR and M of the star-forming galaxies, and to quantify this, 
we computed a 2cr-clipped least squares fit, assuming the rela- 
tion SFR ^ AiM/lO^^Mef. We found values for (A, /?) equal 
to (23.47, 0.70), (38.93, 0.73), (46.26, 0.65), (127.98, 0.85), re- 
spectively in the four redshift bins from z ~ 0.3 to z ~ 2.5. 
The increase in the normalization A is robust, and caused by 
the global increase in the star formation rate with redshift. The 
steepening of the slope, instead, is less robust, since it depends 
on the threshold used for the cr-clipping and on the SED fitting 
parameters, e.g., the minimum galaxy age. With a 3cr-clipping, 
slopes can vary by ~ 20%. Considering these uncertaint ies, our 
results are in broad agree ment with previous studies CDaddi et al.l 
l2007RlEibaFetani2007l) . 

It is interesting to note that the correlation between SFR and 
stellar mass holds also using the SFRfit estimates. While the two 
correlations are consistent at intermediate redshifts, the slope of 
the SFRfit one is steeper at low redshift (0.90) and milder at high 
redshift (0.42), because of the systematic trends shown in Fig. [3] 

4.2. The evolution in the cosmic star formation rate density 

The most concise representation of the evolution in the star for- 
mation rate across cosmic time is the star formation rate density 
(SFRD), which we show in Fig.|5] We show both the IR-derived 
estimates (continuous lines) as well as those obtained purely by 
the SED fit (dashed line). In this figure only, we extend the anal- 
ysis to the redshift bin 2.5 < z < 3. In this redshift regime, the 
IR templates are hardly representative, since the 24 /im band is 
no longer repr e sent ative of the dust emission peak, although the 
iDale & Heloul (1200 2) models are still formally applicable. For 
these reasons, the SFRir+uv estimate is tentative and should be 
interpreted with caution. 

Because of the complex selection criteria that we adopted, 
it was not easy to compute the corrections necessary to include 
the contribution of galaxies fainter than those included in our 
sample. Because of our mid-IR selection, we assumed that we 
include nearly all highly star-forming galaxies. The missed frac- 
tion consists therefore mainly of galaxies with low star formation 
activity, sampled by the z-selection. A good approximation of 
our incompleteness can be obtained by fitting the star formation 
rate function computed on the z-selected sample in each red- 
shift bin. Volumes were computed with the lfV,„ax method. We 
found that the SFR functions can be fitted by Schechter functions 
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Fig. 5. Upper panel: the redshift evolution of the total cosmic 
SFRD. Solid lines refer to IR-based estim ates, dashed lines to 
SEP fitting estimates. The compilation of iHopkins & BeacomI 
( 12006.) is shown as small gray dots. Lower panel: the redshift 
evolution of the cosmic SFRD by galaxies of different masses. 
Only galaxies above the completeness limit in stellar mass at 
each redshift (see last column of Table[T]) are shown. Solid lines 
refer to IR-based estimates, dashed lines to SED fitting esti- 
mates. Black, red and blue symbols correspond to different mass 
ranges as shown in the legend. Error bars include uncertainties 
on the SFR estimates and Poissonian errors. 



of slopes 1.5. We then corrected our observations for these 

contributions at low SFR values by extrapolating from the fitted 
functions. The corrections are small for z < 2.5, of the order of 
10%, consistent with the assumption that the sample is complete 
at high star formation rate values. At z > 2.5, the corrections are 
higher than 50%. 

We first note that the total SFRD (upper panel of Fig. 
|5]) which we derived from our sample (black stars) repro- 
duces well the compilatio n of other surveys completed by 
iHopkins & BeacomI (l2006l) . normalized to a standard Salpeter 
(1955') IMF, and is also in agree ment with Le Floc 'h et alj 
( 200 5), Cap uti et all (l2007l) and.Ro dighiero et al. (2009). 

It is more interesting to study the trend in the SFRD for 
galaxies of different stellar masses (l ower panel of Fig . [S]). 
Previous surveys, such as Combo-17 dZheng et alj 120071 that 
trace d the evolution in the SFRD to z ^ 1) and the GDDS 
one dJuneau et alj|2005i, which first presented the evolution in 
the SFRD in a A'^-selected sample to z ^ 2), demonstrated 
that the contribution to the SFRD by the more massive galax- 
ies (typically those with M > 10"Mq) is negligible at low red- 
shift, and becomes much larger at z >> 1 corresponding to the 
major epoch of formation in the more massive galaxies. Our 
data confirm this picture, with an increase by a factor 20 from 
z - 0.5 to z - 2. The sharp decline at low z suggests that 
these structures must have formed their stars at earlier epochs. 
Keeping in mind the caveat mentioned above about the reliabil- 
ity of SFRiR+uv at z > 2.5, we note that the increase in the star 



Fig. 6. The contribution to the cosmic SFRD by galaxies of given 
mass, in two redshift bins. Solid and dashed lines are computed 
assuming f{M) - 1 . The dashed line in the lower panel corre- 
spond to the z - 1 line of the upper panel. Dotted lines are com- 
puted assuming f(M) + 1 (see text). Black fi lled dots: this work; 
blue sq uares: Combo-17 ([Zheng et al.ll200"7l) : red open triangles: 
GDDS dJuneau et al.ll2005lK 



formation rate density in M > lO^'M© galaxies appears to halt 
at z > 2.5. This is expected since, at these z, the nu mber den- 
sity of M > 10' ^Mq ga laxies drops quickly (iFontana et al. .20061 
iMarchesini e"tani2008h . 

The evolution in the SFRD for galaxies of different stellar 
masses is one of many pieces of evidence of downsizing. Indeed, 
downsizing is evident as an evolution in the slope of the SFRD as 
a function of redshift, which becomes yet steeper as high mass 
galaxies are considered because of the more rapid evolution of 
massive galaxies, wherever a number of physical processes are 
expected to suppress the star formation more efficiently than in 
lower mass galaxies. 

The exact values of the SFRD depend very sensitively on 
the mass ranges adopted. Our mass bin choice is determined by 
the intention of having good statistics in each bin. Keeping in 
mind this warning, our observations seem to support the down- 
sizing picture for the reasons mentioned above. Indeed, in the 
redshift interval 0.3 - 1.5, where all mass samples are complete, 
the SFRD derived from the IR emission increases by a factor 1.8, 
2.75, and 9, respectively for bins of increasing mass. We also de- 
tect a similar steepening in the slope of the relation between the 
average star formation rate per unit mass (which is treated in 
Sect. 14. 3t and redshift between low and high mass galaxies. 

At all observed redshifts, the more massive galaxies do not 
dominate the global SFRD. To some extent, this may follow 
from the exact choice of the mass bins. However, it is easy to pre- 
dict that the major contributors to the global SFRD are galaxies 
immediately below the characteristic mass M*. The number of 
galaxies of given mass N{M) at any redshift can be represented 
as a Schechter function (N(M) oc {M I M*)" exp{-M I M*)). Given 
the tight correlation between stellar mass and star formation rate 
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described above (M oc M^), and assuming that the fraction of 
star-forming galaxies is f{M), the contribution to the cosmic 
SFRD for galaxies in a logarithmic interval dlogM is simply 



d(SFRD(M)) 
dlog(M) 



oc {M/My^'^^^exp{-M/M'-)f{M) 



(3) 



Since a is in the range ^ -1.2/ - 1.4 (iFontana et al.ll2006h 
and /3 - 0.6 - 0.9 (see above, relating to the IR-based es- 
timate), and ignoring for the moment f{M), the shape of the 
d(S FRD(M))/dlog{M) is that of a Schechter function with pos- 
itive slope, which has a peak around M* and decreases at 
log{M) < logiM*). The shape of f{M) is somewhat more un- 
certain. We simply assume that the fraction of active galaxies 
above the characteristic mass is around 0.5 at z ^ 1, around 0.8 
at z - 2, and will probably be higher at lower masses (around 
0.8 at z ^ 1 and around unity at z ^ 2). As a result, it fur- 
ther decreases the distribution of d(S FRD(M))/dlog{M) at high 
masses, without affecting the low mass regime. 

Using the pa rameters of the galaxy mass function obtained 
in iFontana et al.l (|2006), and the SFR-stellar mass correlation 
shown above, we computed the expected d(S FRD{M))/ dlog{M) 
distribution in two reds hift bins, 0.8-1.2 and 1.5-2, where data 
from different surveys (IZheng et all 120071: lJuneau et al.ll2005h 
could be combined. To convert masses and SFRs of Zheng et al] 
(I2007h from a |Chabrier (2003) to a Salpetei (Il 255|) IMF, we 
used a factor of 1.78 (Bundv et al. 2006) and 1.5 (Ferreras et al.' 
l2005h . respectively. Masses and SFRs of Juneau et al.. (.2005|) 
wer e instead r enor malized from a iBaldrv & Glazebrookl (l2003h 
to alSalpeted (W55) IMF using the factors 1.82 (lJuneau et al.l 
l2005h and 2(Hopkins & Beacom 2006). We note that the argu- 
ment presented above is unaffected by the well known discrep- 
ancy between the directly measured mass density and the one 
inferred from the i ntegration of the star formation rate density 
dWiUdns et al.ll200 8). Our point here does not concern the inte- 
grated amount of stars formed, but only the analytical shape of 
the SFRD. 

The result is shown in Fig. |6] where we plot the expected 
d{SFRD{M))ldlog{M)di&mh\xt\onhy assuming both /(M) = 1 
and a f{M) that decreases with both increasing mass and de- 
creasing redshift. Our data a re in good agreement with those pro- 
vided bv lZheng et al.l ( l2007l) . who inferred SFR estimates using 
the s ame method t hat w e adopt to compute SFRir+uv In con- 
trast, |june^ietal] (|2005l) estimated the SFR from the rest-frame 
UV continuum, and obtained values lower than SFRir+uv- The 
differences can probably be attributed to the different selection 
criteria and SFR estimators used in their work. 

This analysis confirms that the major contribution to the star 
formation is from galaxies of masses around or immediately be- 
low the characteristic mass M* . This statement is consistent with 
the evolution in the SFRD for different mass bins that we show 
in Fig. |5] The combination of the steepening in the SFR-M rela- 
tion, the increase in the SFR values with increasing redshift, and 
the lower fraction of active galaxies at low z, is more significant 
than the effect of the decrease in the characteristic stellar mass 
at increasing z. 

4.3. The specific star formation rate and the comparison with 
theoreticai modeis 

In Fig. |7j we plot the relation between the stellar mass and 
the specific star formation rate (SSFR hereafter) for all galax- 
ies divided into redshift bins. To be able to compare our find- 
ings with the Millennium Simulation predictions, we converted 
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Table 1. Average observed IR- and fit-based SSFR at differ- 
ent redshifts for galaxies above the completeness limit in stellar 
mass. 



our masses and SFR to the IChabried (l2003h IMF used by the 
Millennium Simulation. 

We recall that our SSFR estimates are quite low compared 
to the possible estimates made by adopting standard values for 
SFRiR+uv- On the one hand, the conversion between mid-IR 
flux and SFR includes a lowering factor for high SFRs, as de- 
scribed in Sect. 13. H and in Appendix A. On the other side, t he use 
of more recent models of stellar populations (Marastorn l2005l) 
would slightl y decrease the avera ge stellar masses (by 20% on 
average, see iSahmbeni et al.ll2009l) . and hence increase the de- 
rived SSFR. 

First of all, we notice a strong bimodality in the SSFR dis- 
tribution. Two distinct populations, together with some sources 
lying between the two, are detectable, one of young, active and 
blue galaxies (the so-called blue cloud) and the other one con- 
sisting of old, "red and dead", early-type galaxies (red sequence) 
(see also Fig.|4]i. The loci of these two populations are consistent 
with the selection in lSalimbeni et al.l ([2008) between early- and 
late-type galaxies. It is noticeable how these evolved objects are 
akeady in place even at the highest redshi ft. We reserve a m ore 
complete discussion of these objects to Fo ntana et al.l (|2009|) . 

A trend for the specific star formation rate to increase with 
redshift at a given stellar mass is evident: galaxies tend to form 
their stars more actively at higher redshifts. The bulk of active 
sources shifts to higher values of SSFR wi th increasing redshift. 
Our fi ndin gs are in good agreern ent with iPerez-Gonzalez et al.l 
(120051) and lPapovich et al.l (i2006l) . 

A significant fraction of the sample, increasing with redshift, 
is in an active phase. It is natural to compare the SSFR (which 
has units of the inverse of a timescale) with the inverse of the 
age of the Universe at the corresponding redshift f[/(z). We de- 
fine galaxies with M/SFR < tu{z) as "active" in the following, 
since they are experiencing a major episode of star formation, 
potentially building up a substantial fraction of their stellar mass 
in this episode 0. Galaxies selected following this criterium are 
forming stars more actively than in their recent past. 

At 1.5 < z < 2.5, the fraction of active galaxies in the total 
sample is 66%, and their mean SFR is 309 M^yr^K To com- 
pute the total stellar mass produced within this redshift inter- 
val, it is necessary to know the duration of the active phase. 
For this purpose, we use a duty cycle argument and suppose 
that the active fraction of galaxies is indicative of the time in- 
terval spent in an active phase. We adopt the assumption that 
the active fraction is stable within the redshift bin considered. 
The time spanned in the 1.5-2.5 redshift interval corresponds 
to 1.5 Gyr. By multiplying the fraction of active galaxies by 
the time available, we derived an average duration of the active 



^ Indeed, if M = {S FR)i„„xtu(z), where (S FR)i„,, is the star forma- 
tion rate averaged over the whole age of the Universe at the correspond- 
ing z, the requirement SFR/M > {tu(z)y' implies 5 Fi? > {S FR) ,. 
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Fig. 7. Relation between the specific star formation rate and the stellar mass calibrated to a IChabried (l2003h IMF. Larger dots 
correspond to 24 jjm sources with F24^,„ > 83//Jy, while smaller ones show 24 /im detections with lOjjJy < F24ttm ^ 83/iJy. Plus 
refer to 24 jum upper limits. Dotted lines correspond to constant SFRs of 1, 10, 100 and 1000 Mojr The horizontal dashed lines 
indicate the inverse of the age of the Universe at the centre of each redshift interval. Shaded contour levels (at 0.05%, 10%, 50% 
and 80% level) represent the predictions of the lKitzbichler & White! (120071) rendition of the Millennium Simulation. 



phase of 0.99 Gyr. The average amount of stellar mass assem- 
bled within each galaxy during these bursts is measured to be 
the product of the average SFR and the average duration of the 
active phase, and equals to 3.1 x lO^M©, representing a sig- 
nificant fraction of the final stellar mass of the galaxies consid- 
ered. Although quite simplified, this analysis implies that most 
of the stellar mass of massive galaxies is assembled during a 
long-lasting active phase at 1.5 < z < 2.5. It is important to 
remark that this process of intense star formation occurs di- 
rectly within already massive galaxies, and, given its intensity, 
prevails throughout growth episodes due to merging events of 
already formed pro genitors. A similar point is also stated by 
iDaddi et"aLl (l2007bl) . Independent arguments converging on the 
same result are base d on the tightness of the SFR- mass relation 
(Noeske et al.l |2007|) . on the kinematics of disks ( Cresci et alj 
[2009 ) and on the analysis of the accretio n histories of dark mat- 
ter haloes in the Millennium Simulation jGenel et al.ll2008l) . 

To provide a further physical insight in this process, we com- 
pared our results with the predictions of three recent theoretical 
models of galaxy formation and evolution. Our sample is af- 
fected by mass incompleteness, so only galaxies above the com- 
pleteness limit in each redshift bin were considered in the com- 



parison. We note that this limit depends on the redshift bin. In 
Table [1] we report the average SSFR as a function of the red- 
shift bin for the IR- and the fit-based estimates, along with the 
completeness mass cuts computed as described in Fontana et alJ 
(120061) . SSFRs and masses are calibrated to a ISalpeteH (Il955h 
IMF. 

In Fig. |7] we show the pred i ctions of a semi-analytical 
rendition of Kitzbichler & Whit^ (l2007h of the Mil le nnium 
N-body dark matter S imulation (ISprin gel et al.l 120051: 



ILemson & Springell 120061: iDe Lucia & Blaizo^ |2007b, which 
adopts a WMAPl cosmology. We find that the model predicts 
an overall trend that is consistent with our findings. The SSFR 
decreases with stellar mass (at given redshift) and increases 
with redshift (at given stellar mass). In addition, it forecasts the 
existence of quiescent galaxies even at z > 1.5. However, the 
average observed SSFR is systematically under-predicted (at 
least above our mass limit) by a factor ~3-5 by the Millennium 
Simulation. A similar trend for the Mi l lennium Simulation at 
z ~ 2 was already shown bv Daddi et al.l (l2007bl) . 

We als o compared our fi ndings with the semi-analytical 
models of iMenci et al.l (l2006l) (adopti ng a cosmology con sis- 
tent with WMAPl) and MORGANA dMonaco et al.ll2007l up- 
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Fig. 8. Average SSFR in each redshift bin. Black solid lines re- 
fer to IR-based estimates, dashed lines to the SED fitting esti- 
mate. Error bars for the fit-based estimate are larger than those 
of the IR-based one because of the intrinsic parameter degen- 
eracy involved with the SED fitting procedure. In the upper 
panel, we show the comparison with lMenci et al.' ('2006!) semi- 
analytical model (red triangles), which adopts a Salpeter, (Il955h 
IMF. In the lower panel we compare our data with the predic- 
tions of the Millennium Simulation ( Kitzbichler & Whitel[2007h 
(blue squares) a nd MORGANA (green circles), which adopt a 
IChabried 12003') IMF. The comparison is done above the com- 
pleteness limit in stellar mass, which is diff'erent in each redshift 
bin (see text and Table [TJ- 



dated bv lLo Faro et al.l (12009 '). adopting a WMAP3 cosmology). 
They show very similar trends with respect to the Millennium 
Simulation, with only slightly diff'erent normalizations. 

In Fig. |8] we plot the average SSFR as a function of the 
redshift bin for the IR- and the fit-based estimates as well as 
those predicted by the three theoretical models. Once again, we 
only consider galaxies above the completeness limit in mass for 
each redshift bin (see last column in Table[TJ. The trend depicted 
must therefore not be considered an intrinsic trend, since we ob- 
serve different populations at the different redshifts because of 
the mass cuts. Errors in the average SSFR were estimated by a 
Monte Carlo simulation. 

We find that although all models taken into account repro- 
duce the global observed trend, they predict an average star for- 
mation activity that is lower than observed in most of the mass 
regimes. The observed star formation occurring in situ inside 
massive galaxies is higher than predicted by a factor 3-5 for 
the Millennium Simulation an d for MORG ANA, and around 10 
for lMenci et al.l (i2006i) model. iDaddi et all (i2007bl) claimed that 
the star formation rates predicted by the Millennium Simulation 
are as much as one order of magnitude lower than those ob- 
served at z ~ 2. Similar mismatch at these redshifts between 
observed SFRs and those expected by various kind of theoretical 
models, independent of their ph ysica l proce ss implementations, 
were found and discussed by iPavel (l2008l) . To reconcile data 



Fig. 9. Number density of star-forming galaxies as a function 
of redshift and lower SFR. Filled symbols represent our ob- 
servations corresponding to the SFR thresholds shown in the 
legend. Open gray sym bols at z ~ 2 are the predictions of 
lKhochfar& SUM (120081) model. 

and model predictions, iKhochfar & OstrikeJ ([2008') proposed a 
galaxy formation scenario at high redshift that is mainly driven 
by cold accretion flows. Their model allows an increased star 
formation efliciency resulting in a closer agreement with obser- 
vations. 

We then considered the number density of galaxies with 
star formation rates higher than a fixed threshold as a func- 
tio n of redshift, and compared our observations with the model 
of ' Khochfar & Siifl (|2008) at z ~ 2. Figure |9] compares our 
data (filled symbols) with the theoretical predictions (open sym- 
bols). The SFR l imits were chosen to allow comparison with 
iKhochfar & Si S(l2008|)resultfl The agreement is very good for 
objects with SFR > 9QMQyr^\ while the model slightly under- 
predicts the number density of galaxies with higher levels of star 
formation rate. 

A more comprehensive comparison betw een theoretical 
predictions and observations was presented in iFontanot et al.l 
(2009.) . 

5. Summary and discussion 

We have presented a revised version of our GOODS-MUSIC 
photometric catalog for the GOODS-S field. The major new fea- 
ture of this release, on which our scientific discussion is based, 
is the inclusion of 24 jjm data taken from the Spitzer MIPS pub- 
lic images, of which we provide a self-consistent photometry for 
each object in the catalog. We employed a PSF-matching tech - 
nique, performed by ConvPhot software ( De Santis et al.ll2b07h . 
to measure 24 yum photometry by exploiting the high resolution 
z-ACS image used as a prior to detect the objects' positions. This 

Note that they adopt a lChabrieJ ( |2003[) IMF, so the SFR thresholds 
have been renormalized by adopting the conversion factors in Sect. l4.2l 



p. Sandni et al.: Star formation and mass assembly in high redshift galaxies 



13 



allows us to reduce the effects of confusion noise and deblending 
problems caused by the large size of the MIPS PSF. 

We have used this new catalog to study the star-forming 
properties of galaxies up to z ^ 2.5. We have first compared 
the estimates of the SFR obtained from the total IR luminos- 
ity (SFRiR+uv) and from the SED fitting analysis to the overall 
0.3 to 8 i-im photometry (SFRfit), which are the two major es- 
timators of the SFR that have been used so far in high redshift 
galaxies. We found that the two tracers are consistent overall, 
especially in the redshift range 0.3 - 1.5. The overall median ra- 
tio of SFRiR+uv to SFRfit is around unity, a limited number of 
objects appearing to have discrepant results. The agreement be- 
tween the two estimators appears to depend on redshift because 
the scatter increases with redshift. Moreover, SFR^t is slightly 
systematically overestimated with respect to SFRir+uv in the 
lowest redshift bin (0.3 - 0.6). The scatter increases significantly 
at z > 1.5, where the IR-hUV value is systematically larger than 
the one resulting from the SED fitting analysis for bright objects. 
However, these results stem from a systematic trend that holds 
at all redshifts. In galaxies with star formation rates < lOMo/yr, 
the fit-derived SFR is on average higher than SFRir+uv, while 
the opposite holds at SFR > lOOMo/yr. It is at present difficult 
to ascertain the origin of this systematic trend. It can be due to 
systematics that affect the interpretation of the data, either from 
the SED fitting technique as the target galaxies move in redshift, 
or in the templates used in the extrapolation of the mid-IR ob- 
served flux, especially at high redshift, where it samples the PAH 
region. Alternatively, it could have physical origins, such as a 
metallicity trend or a failure in the assumption that a single atten- 
uation law can adequately model the output from a star-forming 
galaxy. 

Keeping in mind these uncertainties in the estimate of the 
star formation rate, we summarize the basic results obtained 
by adopting SFRir+uv, which we assume to be a more reliable 
tracer because it is not affected by dust extinction. 

- We show that, at all redshifts considered here, there is a cor- 
relation between the stellar mass and the star formation rate 
of star-forming galaxies. The logarithmic slope of this cor- 
relation, after applying a cr-clipping algorithm to remove all 
quiescent galaxies, is in the range 0.6 - 0.9, with some indi- 
cations of a steepening with increasing redshift. 

- The SFRD derived from our sample agrees with the global 
trend already depicted by other surveys. When divided ac- 
cording to stellar mass, it shows that more massive galaxies 
enter their active phase at redshifts higher than lower mass 
ones. At z > 2.5, the increase in the SFRD due to the more 
massive galaxies (with M > 10" M©) appears to halt, in 
bro ad agreement with t he expectations of theoretical mod- 
els dMenci et alj l2004t iBower et all l2006l) . This is mainly 
because galaxies of higher mass become extremely rare at 
these redshifts. The increase in the slope of the SFRD with 
redshift of samples of increasing mass seems to support the 
downsizing scenario. 

- At all redshifts, the main contributors to the cosmic SFRD 
are galaxies around, or slightly below, the characteristic stel- 
lar mass M*. 

- Massive galaxies at z ^ 2 are vigorously forming stars, typi- 
cally at a rate of 3OOM0yr"'. A simple duty-cycle argument 
(see Sect. 14.3b suggests that they assemble a significant frac- 
tion of the final stellar mass during this phase. 

- The specific star formation rate of our sample shows a well- 
defined bimodal distribution, with a clear separation between 
actively star-forming and passively evolving galaxies. 



While these results are grossly independent of the particular 
star formation rate estimate, the specific details appear to depend 
on the chosen indicator In particular, the correlation between the 
SFR and the stellar mass still holds using SFRfit, but it is steeper 
at z - 0.5 and it becomes flatter at high redshift, remaining 
similar in the two intermediate redshift bins. As far as the star 
formation rate density and the specific star formation rate are 
concerned, the two techniques used to estimate the SFR provide 
consistent results up to z ~ 1.5. However, significant differences 
are evident at redshift > 1.5, where the IR-based SFRD flattens 
and the fit-based one begins to decline. The average IR-based 
SSFR increases monotonically to the highest observed redshifts, 
while the fit-based one has a turn-over around z ~ 1 .5 and then 
decreases. The trends above directly reflect the correlation be- 
tween SFRiR+uv and SFRfit. 

We used our results for the redshift evolution in the specific 
star formation rate, and its trend with the stellar mass, to inves- 
tigate the predictions of a set of theoretical models of galaxy 
formation in A-CDM scenario. 

On the one hand, these models reproduce the global trend 
that we find in the data - the most important being the increase 
in the specific star formation rate with redshift, and its trend with 
stellar mass. Somewhat surprisingly, however, the average SSFR 
of galaxies in our sample is significantly higher than predicted by 
theoretical models, in most of the mass regimes. Essentially, af- 
ter including a strict completeness limit in stellar mass, we found 
that the typical SSFR of galaxies at a given mass is a factor at 
least ^3-5 higher than predicted by the models we have con- 
sidered here. This mismatch is very clear for massive galaxies 
(M, ^ IO^'Mq) at z - 2 and for less massive galaxies at z - 1, 
where we are able to complete a higher quality statistical analy- 
sis. 

It is not obvious to ascertain the origin of this mismatch. It 
could be caused by a genuine failure of the models. These mod- 
els often quench the star formation in massive galaxies to prevent 
the formation of blue, giant galaxies at low redshift, and to re- 
produce the existence of red, massive galaxies at high redshift. 
The mismatch that we observe could imply that the feedback and 
star formation recipes adopted are too simplified or incorrect. 
Most of the models considered in this work adopt the commonly 
used star formation scenario, where gas uniformly collapsing to- 
wards the galaxy ce ntre forms a stable disk. New processes are 
being explored (e.g.. lDekel et al.| [2009) that involve a more rapid 
formation of galaxies by cold streams and could lead to closer 
agreement with observations. 

Alternatively, the mismatch could be related to the overesti- 
mate of the stellar mass of the typical star-forming galaxy - for 
instance, due to a combination of the overestimate of merging 
events and of the star formation activity in its past history, i.e., 
at redshifts higher than those sampled in this work. A possible 
consequence of this is the overestimate of the stellar mass func- 
tion predicted by the Mi llennium Simulation at high redshifts in 
relation to observations dKitzbichler & Whit3l2007h . 

Additional evidence of an incorr ect treatment of t he sta r for- 
mation proce s ses wa s presented by iFontanot et alj (l2009h and 
ILo Faro et al.l (|2009|) . In their models, these authors found an 
excess of low mass galaxies at z < 2 and faint LBGs at z > 3, 
respectively, which is probably counterbalanced by the suppres- 
sion of the star formation to reproduce the observed evolution in 
the SFRD. 

However, we remark that the interpretation of the observa- 
tions is affected by a number of uncertainties, such as uncer- 
tainties in the stellar mass estimates or uncertainties originating 
from both the templates used to convert 24 /vm fluxes into total 
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infrared luminosities (see Appendix A) and the SED fitting anal- 
ysis. An example of these uncertainties is illustrated by the mis- 
match between the integrated star formation rate density and the 
stellar mass densi ty. This disagreernent could be alleviated by 
an evolving IMF dWilkins et al.ll2008l;lDavell2008h . which would 
provide lower values of SFR. 

To make conclusive, quantitative statements in this direction, 
we ultimately need to improve the reliability of the SFR mea- 
surements, especially for high redshift galaxies. Forthcoming IR 
facilities, such as Herschel and ALMA, will probably provide us 
with a more coherent picture. 
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Appendix A: The estimate of total infrared 
luminosity 

We describe our method to convert mid-IR fluxes into total in- 
frared luminosities and compare the results obtained from differ- 
ent template libraries. 

By assuming that the IR emission is primarily caused by dust 
heating during star formation, the SFR re mains proportion al to 
the dust thermal emission of the galaxy. iKennicuta (Il998l) as- 
serted that the total infrared luminosity Lir emitted between 8 
and 1000 i-im is a good tracer of the SFR. At the redshift of 
our interest (z ~ 0.3 - 2.5), MIPS 24 yum band probes the rest- 
frame mid-IR emissio n, which has been demonstrated to cor- 
relate with Lir (e.g. IChary & Elbaz |l200ll; lElbaz et all 120021; 
iPapovich & BeiilllOOl Takeuchi et al l2005riai7et al.ll200a) . 

Since dust emits most of its light at longer wavelengths, 
Lir extrapolation from 24 fim observations can be a difficult 
procedure. The rest-frame region under study is particularly 
complex given the presence of polycyclic aromatic hydrocarbon 
(PAH) lines. In principle, these aspects could easily lead to non- 
negligible errors in the determination of Lir values. Therefore, 
we employed two different methods to infer Lir, examining and 
comparing both synthetic and empirical templates using our Ks- 
selected subsample (subsample A). 

Firstly, according to p rocedures normall y adopted in the lit- 
erature, we con s ideredlDale & Heloul (12002 *) (DH hereafter) and 
ICharv & ElbazI (1200 lb (CE hereafter) synthetic libraries. These 
models do not extend to sufficiently short wavelengths to fit the 
SED shape, i.e., they do not include the stellar contribution, and 
at z > 1 only MIPS 24 fim band data can be fitted, the other 
bands moving out of the model range or being shifted to a wave- 
length region dominated by star light. Each model is associated 
with a given Lir. In the case of the CE library, each model 
is provided with its absolute normalization, and hence with a 
given total infrared luminosity. As for the DH library, we as- 
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Fig.A.l. Compari son between the total infrared lumino s ity es- 
timated with Dale & Heloul (|2002|) and ICharv & Elba3 (1200 ll) 
templates in different redshift bins. The two synthetic Ubraries 
predictions are in good agreement within a factor < 3. 



colour. In both cases and for each source, once the ^-correction 
had been applied, we selected the model that reproduced the 
24 fim observed luminosity and normalized it using the flux dif- 
ference between the model and the observed galaxy. 
In Fig. lA.ll we show a comparison between Lir predicted by 
the two different model libraries. As already noticed in literature 
dPapovich et alj[2006t iMarcillac et al.l [20061) . they give consis- 
tent results to within a factor 2-3, with the highest differences 
appearing at high redshifts and high luminosities. However, it is 
possible to observe a few trends with both luminosity and red- 
shift that depend on the specific template details. 

A new kind of approach in our work consists of employing 
empirical local spectra to fit the overall galaxy SED shape in- 
stead of the 24 yum luminosity. These spectra include the stellar 
contribu tion as well a s the effects of dust emission. We con- 
sidered PoUetta et all (l2007b spectra (POL hereafter). This li- 
brary comprises early-type, spiral, starburst, and different kinds 
of AGN local templates, ranging from ~ 100 A to ~5 mm. The 
wide wavelength extension allows us to fit the SED shape us- 
ing the multiband catalog. After performing several tests, we de- 
cided to carry out a fit over 5 bands (IRAC bands -i- 24 //m), 
excluding the optical range, where the evolution in the stellar 
component has a more significant effect than the dust emission. 
Since we removed AGN-dominated objects from our sample, we 
only fit early-type and star-forming models. If we instead include 
AGN SEDs, the estimated Lir is only poorly affected, resulting 
in slightly lower values (see also Papovich et al. 2007). Despite 
the poor statistical weight (one amongst five bands) and higher 
noise, the 24 fim band is still fitted acceptably. The mean de- 
viation between observed and fitted magnitudes for 24 fim de- 
tections is (-0.28 + 0.53) mag, 68% of objects being consis- 
tent within 0.48 mag and 95% of objects being consistent within 
1 mag (considering the F24 > 20 //Jy subsample). An exam- 
ple of different kinds of spectral energy distributions fitted by 



16 



P. Santini et al.: Star formation and mass assembly in high redshift galaxies 



0.01 



r ; ' i 

- • ■ ••*. *. : . '-j^, . • \ ■ 

- XX .X '. ' ' ' ' ■ 

0.3 s z < 0.6 


^ 1 1 h 

■ ■ '•■ '?•■'. ■ . ■ • ■ 

"xVxx X • • • 

0.6 s z < 1.0 ; 


j 1 1 1 lll| 1 — 1 1 1 1 ll| 1 — 1 1 1 1 1 

\ 

■ "xx X * ' 

X « 

X " 1.0 S z < 1.5 ; 

- 1 1 1' 


i 1 1 1 llll 1 — 1 1 1 1 1 lll 1 — 1 1 1 

1.5 S z < 2.5 : 

- - 



10 100 

F24[AiJy] 



1000 10 100 

F24[MJy] 



Fig. A.3. C omparison betw een the total in frared lurninosity esti- 
mated with lDale & Heloul f2002) and Poll etta et all ilQQ'm tem- 
plates. Black and gray dots and red crosses indicate objects fitted 
respectively by starb urst, spiral-like and early-type templates by 
iPoUetta et alj ( l2007l) . Synthetic models appear to overestimate 
LiR at bright luminosities in the highest redshift bin when com- 
pared to empirical templates. 



iPoUetta et all ( l2007b templates is reported in Fig. \K2\ 
Figure IA.3I shows the comparison between Lir estimates ob- 
tained with the DH synthetic library and the POL empirical one. 
Source classification based on the spectral shape fitting agrees 
with that assumed from 24 yum emission. As expected, the only 
galaxies fitted by early-type SEDs (red crosses) have very faint 
24 fj.m emission. Moreover, objects tend to be fitted by starburst- 
like models (black dots) rather than spiral-like ones (gray dots) 
as redshift increases. Although the scatter between the two 
Lir estimates is larger than in the previous case (note that, in 
contrast to Fig. lA.ll the y-scale is logarithmic), we found global 
consistency between the two adopted procedures. The only rel- 
evant deviation affects the highest redshift bin: in this redshift 
range synthetic models appear to overpredict, by up to a factor of 
10, Lir for br ight objects. The same behaviour was observed by 
iPapovich et al . (2007), where 24 jum flux, fitted with DH library, 
seems to overestimate the SFR (which is almost proportional 
to Lir) with respect to the case where longer wavelengths (70 
and 160 jim MIPS bands) are considered as well. They corrected 
this trend using an empirical second-order polynomial. A similar 
trend of the 24 jum flux o verestimating L i r at br ight luminosities 
was al so pointed out by iBavouzet et al.l (l2008h andlRieke et all 
(12009^. In our work, we used the iDale & Heloul (|2002|) library 
for consistency with previous works and for which the low- 
ering correction need ed to estimate the SFR was provided by 
IPapovich et all (|2007|) . 

In Fig. IA.4I we show Lir, estimated using both synthetic 
(DH) and empirical (POL) templates, as a function of redshift. 
Once again, the two libraries provide fairly consistent results. 
All in all, the conversion between fluxes and total infrared lu- 
minosity does not depend on the assumed templates in a signif- 



icant way. T he observed LiR-redshift rel ation agrees with that 
presented bvl Perez-Gonzalez et al.l (l2005h . As expected, we are 
only able to detect the more luminous sources as we move to 
higher redshift. 
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Fig.A.4. Relatio n between tot al infr ared luminosity Lir, 
computed using Dal e & Heloul (|2002|) (upper panel) and 
[Polletta et al. (2007) templates (lower panel), and redshift. Black 
dots represent fluxes above 83 /iJy, gray + black dots refer to our 
sample limited at 20 jiiy. 



Appendix B: Error analysis for SFRfu 

In Sect. 13.21 we compared two different tracers for the star for- 
mation rate, one derived from the 24 jim emission and the other 
inferred from a SED fitting technique. Although the two esti- 
mators were found to be consistent overall, in a more careful 
analysis they were found to exhibit a systematic trend depend- 
ing on the SFR itself as discussed in the text. To assess the reli- 
ability of the SFR measurements inferred from the SED fitting, 
and hence concrete evidence of this characteristic trend, we ap- 
plied an error analysis simi l ar to that widely adopted in similar 
cases (IPapovich et £0120011 iFoiTtana et al. 2006. 2009). Briefly, 
the full synthetic Ubrary used to identify the best fit spectrum 
was compared with the observed SED of each galaxy. For each 
spectral model (i.e., for each combination of the free parameters 
age, T, Z, E{B - V)), the probability P of the resulting was 
computed and retained, along with the associated SFR. In the 
case of galaxies with only photometric z, an additional source 
of error was the redshift uncertainty. To account for this, the er- 
ror analysis was completed by allowing the redshift to be a free 
parameter at a local minimum around Zpi,ot- 

We chose two subsamples by performing the error analy- 
sis, for the cases SFRir+uv << SFRjit and SFRir+uv >> SFRfit. 
These two subsamples were called subsample Al and subsam- 
ple A2, and selected to complete a reliable statistical analysis 
in both cases. Subsample Al is made of galaxies having 0.6 < 
z < 1.0, SFRiR+uv < lOMojr-' and SFRiR+uv/SFRfit < 0.8. 
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Fig. B.l. Probability distribution of the star formation rates es- 
timated through the SED fitting technique averaged over each 
subsample. Upper and lower panels show subsample Al and A2 
(see text), respectively. The solid line shows the probability as- 
sociated with the SFR inferred from the generic fitted template, 
SFR,, as a function of the ratio SFR//SFR/,gy/y,7. The dotted line 
shows the distribution of SFRiR+uv/SFRip.,,/,, for all the galaxies 
in the subsample. The dashed line shows the average value for 
SFRiR+uv/SFRi,(,s,/,Y. 



Subsample A2 consists of galaxies in the redshift bin 1.5 - 2.5, 
with SFRiR+uv > lOOM^yr^ and 2 < SFRiR+uv/SFRfit < 20. 
Highly obscured AGN candidates were removed from both sub- 
samples. 

The results of our analysis are shown in Fig. IB. II We 
show the probability distribution of the star formation rates 
estimated by the SED fitting averaged over each subsample. 
We plot the probability associated to the SFR inferred from 
the generic fitted template, SFR,, as a function of the ratio 
SFR,7SFRip.5,/„. We also show the distribution and the average 
value of SFRiR+uv/SFRfe^s,//,. 

As expected, the probability curve is much wider at higher 
redshifts than at lower z, where the spread in our data also be- 
comes larger (see Fig. [3j. This is caused by the faintness of 
the galaxies observed at these redshifts and to the larger k- 
corrections. However, it is clear from Fig. IB. 11 that the SED fit- 
ting for galaxies in both subsamples is well constrained, and 
that the inferred average SFR is significantly higher (lower) 
than the one derived from the IR emission. Indeed, the average 
SFRiR+uv/SFRfegj,/,, (dashed line) as well as their distribution 
(dotted line), occupy the tail of the SFRz/SFRie,,/^,, distribution. 

We also inspected the individual probability curves, consid- 
ering, for each source, the ratio of the probability that the SFR 
from the SED fitting equals SFRir+uv to the best-fit probability 
P(SFR,=SFRiR+uv)/P(SFR,=SFRi,,,/„). The fraction of galax- 
ies with a ratio higher than 0.4 is 17.4% and 7.2% in subsamples 
Al and A2, respectively. 

Finally, according to a Kolmogorov-Smirnov test, the prob- 
ability that the SFR values derived from the two tracers are 




Fig. A.2. Different typical observed spectral energy distributions fitted bv lPolletta et alj (l2007h templates. From top left, in anticlock- 
wise direction, galaxies have been fitted by an elliptical, a spiral and two starbursts templates. The green shaded region indicates the 
wavelength range used for the fitting procedure. 



